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Recently (PRL 96, 106804 (2006)) it was suggested that cobaltocene(CC) molecules encapsulated 
into (7,7) carbon nanotubes (CNT@(7,7)) could be the basis for new spintronic devices. We show 
here based on impact molecular dynamics and DFT calculations that when dynamical aspects are 
explicitly considered the CC encapsulation into CNT@(7,7) does not occur, it is prevented by a 
dynamic barrier mainly due to van der Waals interactions. Our results show that CNT@(13,0) 
having enough axial space for encapsulation but no enough one to allow freely rotation of the 
cobaltocene molecule would be a feasible candidate to such application. 

PACS numbers: 68.65.-k, 61.46.+W, 68.37.Lp, 71.15.-m 



Carbon nanotubes (CNTs) have been one of the most 
studied materials in the last years. Their nanosize and 
diameter have been exploited as a basis for a large va- 
riety of applications [T]. They exhibit very interesting 
electrical and mechanical properties. Among these prop- 
erties, the ability of encapsulating atoms and molecules 
[21 [3] can be used to engineer one or quasi-onedimensional 
systems. 

For instance, the encapsulation of fuUerene molecules 
into CNTs lead to self- assembled structures varying from 
linear chain (generically named peapods) [H |3] to very 
complex ones such as multi-helices [H [51 In gen- 
eral, there is not a strong electronic coupling between the 
fuUerenes and the CNTs due to their weak interactions, 
mainly van der Waals (vdW) ones [7[ [B] • When fuUerenes 
are replaced with structures containing metallic atoms 
more complex structures can be formed with potential 
significant electronic interactions between molecules and 
CNTs [9l[10]. Diameter selective effect for the encapsu- 
lation of cobaltocene molecules into CNT has been ex- 
perimentally reported 

More recently, Garci'a-Suarez and collaborators(GFL) 
[T^ [T51 [H] suggested that systems composed of cobal- 
tocene (CNT@(7,7) or CNT@(8,8)) could be the basis 
for new spintronics devices. Cobaltocenes are molecules 
composed of two aromatic pentagonal rings (C5H5) sand- 
wiching one cobalt atom. The authors based their conclu- 
sions on DFT (density functional theory) calculations us- 
ing the SIESTA [T5J[Tn] code. No temperature or dynam- 
ical effects were considered. Also, it must be considered 
that it is a well-known fact [ITl [18] that DFT methods do 
not describe well the vdW interactions, especially in the 
GGA (generalized gradient approximation) used by GFL 



[12 [12] ■ In general, LDA (local density approximation) 
underestimates the bong-length values while GGA tends 
to overestimate them. In this sense in problems where 
the vdW interactions are very important, as in the case of 
cobaltocene encapsulated into CNTs, LDA would be the 
best choice, because it underestimates the bond- lengths, 
it overestimating the covalent aspects and, consequently, 
captures part of the missing vdW interactions. 

Recent classical molecular dynamics simulations [3 [6] 
of the encapsulation of fuUerenes into CNTs have shown 
that dynamical aspects and vdW interactions are of fun- 
damental importance. In order to determine whether 
this is also the case for the encapsulation of cobaltocenes 
into CNTs we have carried out extensive molecular dy- 
namics simulations and DFT ab initio calculations for 
those systems. We have investigated the CNT@(7,7), 
CNT@(8,8), CNT@(13,0), and CNT@(14,0). We have 
considered cases of frozen and free to relax CNTs. Our 
results show that when the dynamical aspects are explic- 
itly taken into account, the cobaltocene encapsulation 
into CNT@(7,7) does not occur, in disagreement with 
GFL [H results. For the others investigated CNTs 
the encapsulation was observed. Initially we will dis- 
cuss the results from MD simulations. We carried out 
a systematic study of impulse MD calculations, where 
the cobaltocene molecules, at different relative orienta- 
tions and initial velocities (1.0, 5.0, and 10 A/ps) are 
directed toward the CNTs, as schematically shown in 
Figure [T] The MD simulations were carried out using 
the UFF (universal force field), implemented in Cerius2 
package [19]. UFF contains bond stretch, bond angle 
bending, inversion, torsion, rotational, and vdW terms. 
It has been shown to produce accurate results for organ- 




FIG. 1: Different initial configurations for the molecular dy- 
namics simulations. We have considered cobaltocene config- 
urations at parallel (a), 45° tilted (b), and perpendicular (c) 
orientations in relation to the nanotube axial axis. 



ics [201111] and metals |2^. We have successfully used it 
to study carbon nanotubes [22, carbon nanoscroUs [24], 
and organics over metallic surfaces [25] . 

In Figure [2] we present representative snapshots from 
MD simulations for the different nanotubes considered 
here. A better view of the process can be obtained from 
the movies in the supplementary materials |26j . 

For all the cases we investigated (different velocities 
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FIG. 2: Representative snapshots from molecular dynam- 
ics simulations for the CC encapsulation processes of a) 
CNT@(7,7); b) CNT@(8,8); c) CNT@(13,0); d) CNT@(14,0), 
respectively. For all the cases considered here (different veloc- 
ities and initial cobaltocene orientations) we did not observe 
the cobaltocene encapsulation into CNT@(7,7) (a), while for 
all the other nanotubes considered it was observed (b-d). 



and orientations) we did not observe the cobaltocene en- 
capsulation into CNT@(7,7). The cobaltocene molecules, 
even at the most favorable situations, are "trapped" at 
the nanotube borders (see videos 01 e 02). Although the 
configuration of the cobaltocene inside the nanotube is 
energetically favorable, there is a dynamic barrier, mainly 
due to vdW interactions, that prevents the cobaltocene 
encapsulation into CNT(8)(7,7). We have also carried 
non-impulse dynamics placing CC near CNT ends and we 
obtained similar results, no encapsulation was observed 
for CNT@(7,7). 

For the other investigated tubes (CNT@(8,8), 
CNT@(13,0), and CNT@(14,0)) we observed the cobal- 
tocene encapsulation. For the cases CNT(8)(8,8) and 
CNT@(14,0), there is enough free space so the barrier 
for cobaltocene rotation is small. They can move almost 
freely. But, as well pointed out by GFL [T5J [T3] this 
rotation freedom will render the nanodevices useless at 
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FIG. 3: Snapshots from rigid-body molecular mechanics sim- 
ulations of the encapsulation processes of a cobaltocene into 
CNT@(7,7). Even at favorable initial configurations (a-c) the 
molecule never crosses the CNT border (c). 



room temperature, since the potentially magnetic order- 
ing would be destroyed by the available thermal energies. 

On the other hand for the CNT@(13,0), there is enough 
axial space to allow an easy encapsulation, but not 
enough to allow freely rotations (see video 03). Thus, 
from a structural point of view our results suggest that 
CNT@(13,0), and not CNT@(7,7) as proposed by GFL 
[121 El) would be a feasible candidate for these kind of 
applications. 

As mentioned before our conclusions are in clear dis- 
agreement with GFL [HI IS] ones and the origin of these 
discrepancies needs to be addressed. One obvious pos- 
sibility it is that this is just the consequence of differ- 
ent geometrical results, and in this case the reliability of 
our molecular force field results against the DFT-SIESTA 
ones needs to be established. 

In order to do this we carried out a comparative study 
of the geometries of the cobaltocene molecules, the nan- 
otubes, and some selected configurations involving cobal- 
tocene and nanotubes. We contrasted the results from 
UFF with the ones obtained using DMol3 gZl HH] and 
SIESTA TB] code. DMol3 is state of art ab initio DFT 
methodology. In all our DMol3 simulations we have 
used relativistic all-electron DFT total energy approach 
in LDA (Wang-Perdew exchange-correlation func- 
tional) and GGA |30j approximations. A double numeri- 
cal basis set with polarization functions was also consid- 
ered. 

To the Siesta |3TJ calculations we have used the 



TABLE I: Summary of the main geometrical features for the 
investigated CNTs. The lateral (L) and vertical (H) cobal- 
tocene dimensions are also displayed. Results in Angstrons. 



GGA(Dmol3) 
LDA(Dmol3) 
UFF(Cerius) 

GGA(Siesta) 
LDA(Siesta) 



(7,7) 
9.525 
9.489 
9.470 
9.631 
9.540 



10.777 
10.753 
10.740 
10.969 
10.868 



(13,0) 
10.298 
10.218 
10.160 
10.329 
10.226 



(14,0) 
11.024 
11.023 
10.934 
11.111 
11.009 



L 
4.602 
4.596 
4.535 
4.442 
4.452 



H 
3.486 
3.360 
3.124 
3.552 
3.443 



TABLE II: Difference (in Angstroms) between the different 
methods for the magnitudes indicated in Table I. 

(7,7) (8,8) (13,0) (14,0) L H 

GGA/LDA(Dmol3) 0.036 0.024 0.080 0.001 0.006 0.126 

GGA(Dmol3)/UFF 0.055 0.037 0.138 0.090 0.067 0.362 

LDA(Dmol3)/UFF 0.019 0.013 0.058 0.089 0.061 0.236 

GGA(Siesta)/UFF 0.019 0.013 0.058 0.089 -0.093 0.428 

LDA(Siesta)/UFF 0.070 0.128 0.066 0.075 -0.083 0.319 

GGA(Dmol3/Siesta) -0.106 -0.192 -0.031 -0.087 0.160 -0.066 

LDA(Dmol3/Siesta) -0.051 -0.115 -0.008 0.014 0.144 -0.083 



standard double zeta plus polarization(DZP) basis with 
an energy shift of 0.27 eV to represent the pseudoatomic 
confinement for all atoms. A cutoff of 180 Ry for the grid 
integration was utilized to represent the electronic charge 
density for both (LDA - CA) [33] and generalized gradi- 
ent approximation (GGA - PBE) [34] for CC molecule 
and CNTs. All atoms position were set free and the 
structure was relaxed for forces below to 0.04 eV/A us- 
ing a conjugate gradient algorithm. The pseudopotential 
was constructed according to Troulier Martins scheme. 

In Tables I and II we present a summary of the main 
geometrical features obtained with the different meth- 
ods. As we can see there is an excellent agreement be- 
tween the geometrical data obtained with the classical 
force field and the ab initio DFT ones. In particular, al- 
though the cobaltocene's H(high) dimension using UFF 
is a little off in relation to the DMol3 ones the relevant 
magnitude determining the encapsulation is L, and in this 
case it is in very good agreement with DMol3 results. It 
should be stressed from the experimental point of view 
it is difficult to determine the chirality of CNT based 
on the diameter estimations. Experimental diameter es- 
timatives from transmission electron microscopy(TEM) 
measurements contain errors as large as 30% for CNT 
with diameters less than 1.0 nm and about 10% for CNT 
of large diameters [35] . Considering these error bars it is 
not possible to discriminate (7,7) from (13,0) CNTs. 

In Table III we present the differences of the "free 
space" (not considering the excluded volume due to the 
van der Waals repulsions) for the different methods. Al- 
though the differences are very small it could be argued 
that if we are in the limit cases these differences could be 
decisive determining whether a molecule would be encap- 
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LDA(Dmol3) 
GGA(Dmol3) 
UFF(Cerius) 
LDA(Siesta) 
GGA(Siesta) 



(7,7) 

4.895 
4.923 (0.028) 
4.935 (0.040) 
5.088 (0.193) 
5.189 (0.294) 



(8,8) 

6.159 
6.175 (0.016) 
6.205 (0.046) 
6.416 (0.257) 
6.527 (0.368) 



(13,0) 

5.624 
5.696 (0.072) 
5.625 (-0.082) 
5.774 (0.150) 
5.887 (0.263) 



(14,0) 

6.429 
6.422 (-0.007) 
6.399 (-0.030) 
6.557 (0.128) 
6.669 (0.240) 



TABLE III: Difference (in Angstroms) between the tube diameter value and the lateral (L) cobaltocene dimension, for the 
different methods considered here. In parenthesis are displayed the differences relative to the LDA ones. 



sulated or not. In order to rule out this possibility we run 
a final test where we used the SIESTA-GGA geometries 
for the tube and cobaltocene molecules in a simulation 
of rigid body (the geometries are keep fixed during the 
process of energy minimization) . Even in this limit case 
using the SIESTA geometries the cobaltocene molecule 
is prevented to be encapsulated by the vdW interactions 
(not well described in DFT approaches) (Figure [s]) . It 
should stressed that for CNT diameters considered here 
no significant charge transfer was observed from SIESTA 
calculations These results strongly indicate that the 
encapsulation of cobaltocene into CNT@(7,7) is not pos- 
sible, in contrast with GFL results [12 . 

In summary, we have investigate using impact classi- 
cal molecular dynamics and ab initio DFT methods the 
encapsulation processes of cobaltocene molecules (CC) 
into carbon nanotubes. In contrast with previous DFT 
studies [HI [T3] our results show that, the CC encapsu- 
lation into CNT@(7,7) it is not possible when dynamical 
effects are explicity taken into account. Our results also 
show that CNT@(13,0) would be a better candidate for 
spintronic applications. These results also point out that 
conclusions based on DFT results of geometrical config- 
urations where van der Waals are of central importance 
should be take with caution. 
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